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qq ' Abstract 

O" 

We consider the N = 4 SYM theory in flat 3 + 1 dimensional spacetime with a time 

dependent coupling constant which vanishes at t = 0, like g\ M = t p . In an analogous quantum 

mechanics toy model we find that the response is singular. The energy diverges at t = 0, for 

a generic state. In addition, if p > 1 the phase of the wave function has a wildly oscillating 

behavior, which does not allow it to be continued past t = 0. A similar effect would make 

the gauge theory singular as well, though nontrivial effects of renormalization could tame this 

singularity and allow a smooth continuation beyond t = 0. The gravity dual in some cases 

is known to be a time dependent cosmology which exhibits a space-like singularity at t = 0. 

Our results, if applicable in the gauge theory for the case of the vanishing coupling, imply that 

the singularity is a genuine sickness and does not admit a meaningful continuation. When the 

coupling remains non-zero and becomes small at t = 0, the curvature in the bulk becomes of 

order the string scale. The gauge theory now admits a time evolution beyond this point. In 

this case, a finite amount of energy is produced which possibly thermalizes and leads to a black 

hole in the bulk. 
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1 Introduction and Summary 

The resolution of singularities is an outstanding problem in the study of gravity. The gauge 
theory/gravity correspondence provides a non-perturbative framework for the study of gravity, 
and one would hope that it can shed some light on this question. 

With this motivation some cosmological solutions which admit natural gauge theory duals 
have been constructed and studied recently [1, 2, 3, 5, 6]. These solutions can be thought of as 
deformations of AdS^, x S* 5 (in Poincare coordinates) and are dual to strongly coupled M = 4 
Yang Mills theory in flat 3 + 1 dimensional space-time subjected to external time dependent or 
null sources. For other classes of solutions of this type, see [7], [8]. A related approach to gauge 
theory duals of cosmological singularities has been pursued in [10] and more recently in [11]. 
Ideas about finding signatures of spacelike singularities inside black holes in the dual gauge 
theory are described in [9]. 

In this paper our interest is in the time dependent cases. These have been further studied in 
[3]. Here, the only source is a time dependent coupling of the gauge theory. At early times, the 
'tHooft coupling in the gauge theory is large and varies slowly. The AdS/CFT correspondence 
tells us that the gravity dual is a non-normalizable deformation of AdS§ x S 5 sourced by the 
dilaton field. As time evolves, the gauge theory state evolves in response to the time dependent 
coupling. On the gravity side, the background evolves according to the supergravity equation 
of motion, subject to appropriate boundary conditions. In particular the dilaton starts back 
reacting leading to a nontrivial metric. 

Our main interest will be in situations in which the dilaton, e*, vanishes at time t — 0. In 
the corresponding gravity solution a space-like singularity appears at t — which extends all 
the way to the boundary 2 . One would like to know if this singularity is a genuine sickness in 
the theory or if it merely signals a breakdown of the supergravity approximation. Since the bulk 
theory has a boundary dual which is formulated in a precise fashion, we can ask this question 
in the dual theory. One would like to know if the boundary theory is sick at the singularity or 
if it allows for a continuation past the point where the dilaton vanishes 3 . 

In this paper we try to analyze this question in some detail. Prompted by the cosmological 
solutions we ask the following general question first: Consider the M = 4 Yang Mills theory 
subjected to an external time dependent dilaton. We take the dilaton to be of the form, 

e* = (~ty, t < 0, (1) 



2 However, the 4 dim. metric as seen in the gauge theory is flat, after a conformal transformation, as discussed 
in section 5. 

3 More precisely, the supergravity approximation breaks down before the singularity forms, once the curvature 
gets to be of order the string scale. Using the boundary theory we would like to find out if there is a continuation 
past this region. 



so that it vanishes as t — > 0~. In eq.(l), the index p can take any positive real value. In such 
a situation the question we ask is whether the response of the gauge theory to this external 
source is singular or not, as t — »■ 0~ ? 

It is useful to first consider a quantum mechanics model which has an analogous coupling 
to the dilaton. We find that the response in this quantum mechanical system is singular, for 
all values of p > 0, in the following sense: The time dependent dilaton pumps energy into the 
system, and as t — »■ 0~ we find that the energy diverges. The nature of the singularity does 
depend on the index p. For p > 1, when the variation of the dilaton is more rapid, the wave 
function of the system acquires a time dependent phase factor which becomes wildly oscillating 
and diverges as t — ► 0~. As a result the wave function of the system (in the Schrodinger picture) 
does not have a well defined limit as t — »■ CT . In contrast for p < 1 the phase factor does not 
diverge and the wave function has a well defined limit as t — ► 0~. Even so, the energy diverges 
as t — ► CT (this also happens when p > 1). 

The result that for p > 1 the wave function becomes singular, without a well defined limit, 
holds regardless of the state of the system. On the other hand, the conclusion that the energy 
diverges is true for a generic state. For special states in which an appropriate matrix element 
vanishes, the energy can remain finite, as we discuss below 4 . 

The analysis in the quantum mechanics model is most conveniently carried out in the 
Schrodinger picture. The conclusions stated above follow from the fact that near t — the 
potential energy term in the Schrodinger equation dominates the time evolution and the Kinetic 
energy term is sub dominant. 

One can carry out a similar analysis in the gauge theory. Once again if we assume that the 
potential energy dominates near t — one finds that the behaviour is completely analogous 
to that in the quantum mechanics model disussed above. The energy diverges as t — ► CT and 
for p > 1 the wave function acquires a wildly oscillating time dependent phase which does not 
have a well defined limit as t — *■ CT. 

However, in the field theory case, we have not been able to establish that the approxima- 
tion leading to the conclusions above definitely holds. Differences with the quantum mechanics 
model arise due to the infinite number of modes in field theory. These have to be dealt with 
carefully by regulating the theory after introducing a cutoff and incorporating the effects of 
renormalisation in this cutoff theory. We have not carried out this procedure adequately to de- 
termine whether the approximation mentioned above of the potential energy dominating holds, 
and our conclusions about the gauge theory response being singular are therefore tentative 
and not definite. As we discuss below, higher loop effects could sum up to render the kinetic 
term important and tame the singularity, resulting in finite energy production and a smooth 



fluctuations in the energy will diverge in these special states, even if the expectation value of the energy 
remains finite. However such fluctuations are suppressed at leading order in N. 



continuation of the wavefunctional beyond t — 0. Studying this further will require both a bet- 
ter understanding of the calculational aspects of renormalisation in the time dependent gauge 
theories at hand, and a better understanding of the conceptual issues involved in incorporating 
these effects of renormalisation in the Schrodinger picture which we use in this paper 5 . 

The AdS cosmologies described in [1, 3] correspond to the value p = a/3 in eq.(l). From the 
discussion above it follows that, if our approximation of a dominant potential energy continues 
to hold in the gauge theory, the resulting singularity is a genuine sickness which does not admit 
a well defined continuation. Since p > 1 in this case, the wave function in the gauge theory 
description becomes wildly oscillating without a well defined limit as t — ► 0~. Also the energy 
should diverge as t — ► 0~. 

Our analysis shows that the singularity in the bulk arises due to two related reasons. First, 
the dilaton vanishes, resulting in the string frame curvature blowing up. Second, an infinite 
amount of energy is dumped into the system by the vanishing dilaton, resulting in a singular 
back reaction. 

We do not know the explicit form of the bulk solutions whose boundaries are conformally 
flat for values of p, other than a/3- However such solutions should exist since we are specifying 
the dilaton field on the boundary for all times. It is worth pointing out that the singular 
behavior we find is not tied per se to the non-analyticity of the dilaton as t — ► 0, and occurs for 
all values of p > 0, integer and non-integer. Rather, the singular behavior is related to the rate 
at which the dilaton vanishes. As was mentioned above, with our approximations, the behavior 
is more singular for p > 1, when the dilaton vanishes more rapidly, than it is for p < 1. 

A few more comments about the analysis in the gauge theory are also worth making. At 
first sight one might think that when the dilaton vanishes the gauge theory becomes weakly 
coupled and can be analyzed in perturbation theory. This turns out not to be true. Starting 
from a generic state, one finds that the time dependent dilaton excites the fields to large enough 
values so that the cubic and quartic interaction terms are non-negligible near the point where 
the dilaton vanishes and perturbation theory is not valid. This is the essential reason why the 
analysis gets complicated. Based on the quantum mechanics model we find that the Schrodinger 
picture is particularly useful in analyzing the resulting behavior. A WKB-like approximation 
can be formulated in this picture in the vicinity of the vanishing dilaton. This allows the 
leading behavior near the singularity to be analyzed without having to resort to perturbation 
theory. As was mentioned above, in the gauge theory, we have not been able to establish that 
this approximation is indeed correct and thus our conclusions should be taken to be indicative 
rather than definitive. 

It is worth emphasizing that if the potential energy term continues to dominate near the 



3 We thank David Gross and the referee for emphasising these points to us. 



singularity, our conclusions about the gauge theory in the presence of a coupling which truly 
vanishes are valid at finite N and finite 9ym^ an d are therefore a result about the dual closed 
string theory in the presence of stringy and/or quantum corrections. As we will see, the 
behavior of the wavefunctional near the time of vanishing coupling is essentially determined by 
the coupling constant and not the details of the lagrangian. 

It follows from our analysis that to find cosmological solutions which are not sick the dilaton's 
behavior at the boundary has to be modified so that it does not vanish. Once this is done, for a 
smoothly varying dilaton, one does not expect the gauge theory to be singular 6 . In a situation 
where the dilaton profile is chosen to be constant in the far future, reaching a value such that 
the 'tHooft coupling in the boundary theory is large, one expects that the supergravity solution 
becomes a black hole in the far future. This is based on the expectation that the dual gauge 
theory will generically have some non-vanishing energy density in the far future and this energy 
will eventually thermalize. It is worth noting that if this expectation is met, the spacetime which 
is highly curved when the dilaton is small, eventually evolves to a smoothly varying spacetime 
outside the black hole horizon. The fate of the theory will be similar if renormalization effects 
tame the gauge theory even in the case where the coupling truly vanishes. 

An important question which we cannot address here is how the thermalization process 
depends on the dilaton time profile. To answer this question one needs a better understanding 
of the system when the 'tHooft coupling is of order one, for example. 

One might wonder if the formation of the black hole can be avoided in a situation where 
the dilaton profile is time reversal invariant. In such a case states should exist which evolve in 
a time reversal invariant manner. Classically, in such a state all velocities have to vanish at 
t — 0; quantum mechanically, the wave function has to be real at t — 0. Starting from t — 
in such a state and evolving into the future one expects that a black hole will typically form, if 
some net energy is input into the system. Thus we do not expect the absence of a black hole 
in such states, rather these states will correspond to starting with a black hole in the far past 
and ending with one in the far future. In some situations, this conclusion might be avoided, 
but we do not understand at the moment how to identify them 7 . 

The breakdown of perturbation theory near t — 0, in the time dependent case, is quite 
different from what happens in the null-dependent solutions which were studied in [1, 2, 5]. In 
the null case, the effects we describe in this paper are absent, there is no particle production 
and perturbation theory is possibly applicable for correlators of fields at light front times when 



6 We are assuming here that the gauge theory does not have any phase transitions as the dilaton is varied. 
If this assumption is wrong the response of the gauge theory to a smoothly varying dilaton need not be non- 
singular. Such a phase transition can be avoided by working on S 3 , at finite N. We thank S. Wadia for 
emphasizing this point to us. 

7 The case where the boundary theory is on S 3 rather than i? 3 is more promising in this respect, since the 
formation of a black hole then requires the temperature to be bigger than that of the Hawking Page transition. 



the coupling is small. 

We discussed cosmological solutions above with p = \/3, which were studied in [1, 3]. These 
correspond to symmetric Kasner-like solutions, where all spatial directions expand at the same 
rate, or to FRW solutions. There are other asymmetric Kasner solutions which were also 
constructed in [1, 3]. The dual gauge theory in these cases does not live in flat space but 
instead in a space-time with a curvature singularity, which occurs at the time coincident with 
the bulk singularity. Since the boundary metric is non-dynamical, it is difficult to see how time 
evolution on the boundary can be continued past the singularity in such circumstances leading 
to the conclusions that the boundary hologram is sick in these cases as well. 

The behavior of the cosmological solutions near the singularity has some degree of univer- 
sality. For example in the symmetric Kasner and FRW solutions, mentioned above, the spatial 
curvature is different but this feature is irrelevant near the singularity where the evolution of the 
spacetime is determined by the diverging dilaton stress tensor. One expects this to be a more 
general feature- some differences among solutions should become irrelevant leading to the same 
behavior near the singularity. Our conclusions obtained from studying the gauge theory dual 
to the cosmological solutions, within our approximations, do not require any particular state, 
and apply quite generally. One therefore expects these conclusions to hold for all solutions with 
the same common behavior near the singularity, e.g., for all conformally flat four dimensional 
metrics as discussed in section 6. 

To explore this further we discuss Bianchi IX type cosmologies in the presence of the dilaton 
towards the end of this paper. These solutions are constructed by taking solutions of 4 dim. 
gravity coupled to dilaton and embedding them in 5 dim. AdS space as discussed in [1, 3]. 
Interestingly, once the dilaton is excited only a finite number of BKL oscillations occur between 
Kasner regimes. With each oscillation the importance of the dilaton stress energy grows, at 
the expense of the spatial curvature. Eventually all solutions enter a Kasner regime, where 
all directions shrink, and the spatial curvature is relatively unimportant. In a whole family of 
solutions this final Kasner regime corresponds to the symmetric Kasner solution with p = a/3, 
mentioned above. For other solutions the final Kasner regime is one where the three spatial 
directions shrink in an asymmetric manner. The holographic dual for these cosmologies can be 
constructed. Our analysis of the singularity in the symmetric and asymmetric Kasner solutions 
then applies to all these solutions. It will be interesting to explore these solutions and their 
holographic duals further. 

This paper is organized as follows. In section 2 we discuss the gauge theory and introduce 
a toy field theory which captures the essential features of the gauge theory in the dilaton back- 
ground. We also introduce a quantum mechanics model which has many features in common 
with these field theories. The quantum mechanics model is then analyzed in considerable detail 
in section 3. The implications for the toy field theory and the gauge theory are discussed in 



section 4. The connections to the cosmological solutions are discussed in section 5. Other 
Kasner and BKL-like solutions are discussed in section 6. 

Several important details are in the appendices. Appendix A discusses the coupling of 
the dilaton to the Yang Mills theory. Appendix B discusses the time dependent Harmonic 
oscillator. Subleading corrections to the energy, in the vicinity of the vanishing dilaton, are 
discussed in Appendix C. Particle production in the quadratic approximation, for a modified 
non- vanishing dilaton profile, is discussed in Appendix D. The behavior of the Yang Mills theory 
in the presence of a dilaton which varies with Milne time is studied in Appendix E. Finally, 
some discussion about the universal behavior near singularities and about BKL cosmologies in 
the presence of the dilaton, is contained in Appendix F. 

2 The Gauge Theory and a Toy Model 

We will consider the N = 4 gauge theory defined on a flat 3 + 1 dimensional space-time, which 
is regarded as the Poincare patch boundary of an AdS cosmology - e.g. the ones discussed in 
detail in [3] . The bulk dilaton is equal to the coupling constant of the Yang Mills theory, 

e m = 9 2 YM . (2) 

Thus the varying dilaton gives rise to a varying Yang Mills coupling. 

After suitable field redefinitions the Lagrangian for the M = 4 theory takes the form, 

L = Tr{-^F, U F^ + 1 -(D,X a y - 1 -e^[X a ,X b }) 2 

-wr»D^ + e* /2 yr a [x a ,y}}. ' (3) 

There are six scalars, X a , a = 1, • • • 6. And 4 two-component Weyl fermions of 5*0(1, 3), which 
have been grouped together as one Majorana Weyl Fermion of S0(1, 9). The Gamma matrices 
r M , n — 0, 1, • • • 3, and T a , a = 1, • • • 6, together form the 10 Gamma matrices of SO(l, 9). The 
scalars and fermions transform as the adjoint of SU(N). The covariant derivative of the scalars 
is, 

D„X a = d lt -i[A lt ,X a ], (4) 

and similarly for the fermionic fields. The field strength is, 

F^ = d»A v - d v A^ + i[A lx , A v \. (5) 

In appendix A we discuss the Lagrangian above, especially the dilaton couplings, in more 
detail. Let us make one comment here. It is sometimes stated that the dilaton couples to the 
Lagrangian of the M = 4 theory, [19]. The more correct statement is that the dilaton couples 



to the on-shell Lagrangian [20]. This differs from the operator obtained from eq.(3) but only 
by a total derivative term involving the scalars. 

At the singularity, e* goes to zero. As a result, the prefactor in the gauge kinetic energy 
term in eq.(3) blows up. This is the essential complication which must be dealt with in our 
analysis of the time dependent situation at hand. To study it further it is useful to introduce 
a toy model consisting of a single scalar field X with Lagrangian 

L = -±l 1 -(dXf + X% (6) 

Note that the quartic term is "right side up", so that for a constant dilaton this model 
would have a stable minimum. We will see that this model captures all the essential features 
of the real problem of interest. 

In fact it is useful to simplify the model further and to consider a quantum mechanics 
system, with action, 

S = J dt^X - l -u 2 X 2 - X^} (7) 

Starting from the field theory eq.(6) such a system would arise if we keep only one fourier mode 
of momentum k with u 2 = k 2 , and only the quartic self interactions of this mode. 

It will be useful in the discussion below to first consider the quantum mechanics system and 
then return to the field theories eq.(6) and eq.(3). 

3 Analysis of Quantum Mechanics Model 

We now turn to analyzing the quantum mechanics model eq.(7) further. 

It is useful to begin by carrying out a field redefinition which gives rise to a variable with 
a canonical kinetic energy term and analyzing the system in terms of this new variable. We 
define the variable X to be 

X = e-* /2 X. (8) 

Upto a surface term the Lagrangian is 

L = \x 2 - l r \t)X 2 - e*X\ (9) 

Here oj 2 {t) is a time dependent angular frequency that arises due to the time dependent dilaton 
and is given by 

uj 2 (t) = -i\(i>) 2 -h] + ul (io) 

where the dot superscript indicates a time derivative of the dilaton. For a dilaton dependence 

e* = (-i)V<0, (11) 



we have 



with, 



u 2 (t) = -^ + u; 2 (12) 

« 2 = |(p + 2). (13) 

Note that for sufficiently small time, u 2 becomes negative, and the variable X has a tachyonic 
mass term (negative (mass) 2 ), which arises due to the time dependent dilaton. In fact uu 2 
diverges as t — *■ 0, this will be important in the subsequent analysis. 

Let us briefly outline the detailed analysis that follows. Since the X 4 term in eq.(9) is mul- 
tiplied by e*, one might at first expect that this term is not important near t — 0. Accordingly, 
we first neglect it and analyze the resulting quadratic theory. It turns out that due to the 
diverging tachyonic mass the system quite generally gets driven to \X\ — ► oo as t — > 0. As a 
result, the quantum mechanical description in terms of the X variable is not complete. If we 
want to know what happens as t — > 0, and beyond, one needs additional information about the 
behavior at X — ► oo. The diverging value of X also means that the quartic term cannot be 
neglected. 

At this stage it is worth remembering that the field redefinition, eq.(8), is singular at t — 0. 
In fact, we find that the rate at which X diverges is exactly balanced by the rate at which the 
dilaton vanishes, leaving X to be finite, as t — ► 0. This motivates us to study the system in 
terms of the X variable. Although, as mentioned above, the analysis in terms of the X variable 
has already revealed that the quartic term cannot be neglected, we ignore it at first, to gain 
some understanding of the system. Our analysis shows that the system is singular in a manner 
we have described in the section 1. We then incorporate the quartic terms and find that all the 
essential conclusions about the singular nature of the response go through unchanged. 

3.1 The X Description 



As was mentioned above, to begin we drop the quartic term in eq.(9). this gives rise to the 
quadratic action, 

S = j dt]-[X 2 -uj{t) 2 X 2 } (14) 

The quadratic theory can then be analyzed in standard fashion by expanding the field, X, in 
terms of normal modes, 

From eq.(14) it follows that the corresponding operator, X, in the Heisenberg picture, 
satisfies the equation 

X + u 2 {t)X = 0. (15) 

Let us define fit) to be, 

/(*) = ^fV^(-cuot), (16) 

9 



where, 

and H^(x) is the Hankel function which asymptotically, asx^ oo, behaves like 

H W( X ) - t /A e W-(-+|)f)]. (18) 



Then it follows from the standard properties of Hankel functions that /(£) satisfies the 
equation, 

/(*)+w 2 (t)/(*)=0, (19) 

with boundary condition, 

/(t)^e- iuJ(,t ,t^-oo. (20) 

The solution to eq.(15) now is, 

X = -±=[af(t) + (ayr(t)]. (21) 

The momentum conjugate to X is 

P = X = -l=[af(t) + (ayr(t)]. (22) 

The operators, X, and P, satisfy the canonical commutation relation iff, a,^, satisfy the 
standard relation, 

[0,0^ = 1. (23) 

Classical solutions to eq.(15) take the form, 



X(t) = j(-t)[AJ v (-t) + BN„(-t)\, (24) 

where J u , N v , stand for the Bessel and Neumann functions. The constants A, B, are determined 
by the initial conditions, at t — ► —00. For generic initial conditions, B 7^ 0, and as t — *■ _ , 

X(t) ~ (-t)5-^ = (-t) ( ^ } -► 00. (25) 

Here we have used eq.(17), and the fact that p > 0. Thus we see that due to the negative and 
diverging value of u 2 near the singularity, a generic trajectory gets driven out to infinite values 
of the position coordinate. 

Classical states correspond to coherent states in the quantum theory. We see that the center 
of the wave packet for a generic coherent state runs away to infinity due to the tachyonic mass 
term. This shows that further data is needed to make the quantum theory in terms of the 

10 



variable X, and Lagrangian, eq.(14), well defined. This additional data should specify what 
happens to the wave packet once it gets to large values of the X coordinate. 

It is also illuminating to calculate the wave function of the ground state in terms of the X 
description. As t — *■ — oo, the system becomes a conventional Harmonic oscillator with constant 
angular frequency u . Consider the state specified by the condition, 

o|0 >= 0, (26) 

which we will refer to as the ground state. We are interested in asking how expectation values in 
this state evolve with time. We have been working in the Heisenberg picture above. It is useful 
to answer this question by constructing the wave function in the Schrodinger picture. The 
resulting time dependence of the wave function carries information about the time dependent 
expectation values for all operators in this state. 

As discussed in appendix B the Schrodinger picture wave function for the ground state is 
given by, 

^( X) t) = ^^M ) *^\ (27) 

where A is a time independent constant which is fixed by requiring that the state has unit norm. 
Two features of the resulting behavior of this wave function near t — are worth commenting 
upon. 

First, as discussed in Appendix B, it follows that the probability density \ip(x, t)\ 2 is given 

by ' 

|^,0| 2 = ^!e- [ T^ ] . (28) 

It follows from eq.(16), eq.(17), and the properties of Hankel functions that 

\f\ 2 ~(-tyv^oo,t^0-. (29) 

Thus the probability distribution in X become infinity spread out, as t — > 0. We saw above that 
for generic coherent states the center of the wave packet runs off to infinite X. The vacuum is 
a non-generic coherent state, for which this does not happen. The expectation value of < X > 
vanishes in this state, this corresponds to A, B in eq.(24) both vanishing. However we see now 
that even for this state the spreading of the wave function, which is a quantum effect, makes 
the wave function sensitive to large X. 

Second, the exponential factor in eq.(27) gives rise to a phase factor, 

e '[tfr4] „«,[*«£]. (so) 

This phase factor oscillates "wildly" near t — 0. As a result ip(x, t) does not have a well defined 
limit as t — ► 0~. This feature will be crucial in the subsequent analysis that follows. 

11 



We have neglected the quartic interaction term of eq. (9) in the analysis above. We now see 
from eq.(25) that X diverges, as £~ p / 2 , and thus e®X 4 goes like, t p t~ 2p — ► oo, as t — ► 0, and also 
diverges. This means that despite the vanishing dilaton the quartic term cannot be regarded 
as a small perturbation. 

It is worth recalling now that the theory we started with was formulated in terms of the 
X variable, eq.(6). The change of variables from X to X is in fact singular at t — where 
the dilaton vanishes. Moreover, as we will see below, X is in fact finite, as t — > 0, since the 
the rate at which X diverges in eq.(8), is exactly the rate at which e -$ / 2 also blows up. It is 
therefore worth constructing a description directly in terms of the X variable. The resulting 
behavior at finite X will also provide information about what happens at infinite X, which is 
the additional data we seek. We turn to this next. 

3.2 The X Description 

As was mentioned above, to gain some understanding in the X variable we first begin by 
neglecting the quartic term. The action in terms of this variable then takes the form, 

S = J dt e-*]-[X - ulX 2 }. (31) 



Classical solutions take the form 






X(t) = e~X{t) = e*y/(-t)[AJ„(-t) + BN„(-t)]. (32) 

It is easy to check from eq.(17) and eq.(ll) that as t — > 0, 

X(t) ~ Be^N v (-t) -► constant. (33) 

Thus classical trajectories do not reach \X\ — ► oo, but instead are at finite values as t — ► 0. 

From the relation, eq.(8), and eq.(21), it follows that the operator X in the quantum theory 
is 

X = -j^[af(t) + tfnt)]. (34) 

The ground state satisfies the condition given in 8 eq.(27). The resulting wave function in the 
Schrodinger picture (see appendix B) is, 

if,(x, t) = . A e^lr+tl^}. (35) 



8 Note that for the dilaton dependence given in eq.(89), the Lagrangian in terms of the X variable does not 
reduce to that of a standard harmonic oscillator as t — > — oo. This is in contrast to the solution, cq.(84), where 
we do get a standard harmonic oscillator, as t — > — oo, as discussed in appendix E. It is this latter case that is 
better defined in any case, as was discussed in section 3. The ground state in this latter case, which becomes 
the vacuum of the standard harmonic oscillator in the far past, behaves similarly to the vacuum state consider 
here near the singularity. 

12 



The probability to find the system between x,x + dx, is \ip(x, t)\ 2 , and is given by, 



\A 



2 



(x : t)\ 2 = -^^e W^ J . (36) 



From eq.(29), eq.(ll), it follows that |/| 2 e* goes to a constant as t — > 0. Thus \ip(x, t)\ 2 becomes 
a well-defined smooth Gaussian function in the limit t — ► 0. The absolute value of the wave 
function, \ip(x, t)\ thus has a smooth limit as t — ► 0~. Contrast this with the phase of the wave 
function. As discussed in appendix B, 

e -»[(/)*+|]->i/(-tr\ (37) 

as t — > 0~ , and thus the phase of the wave function goes like, 

e ^)*+f]^ _> e i[ (^H (38) 



where C is a constant. Note that for p > 1 the phase factor diverges 9 . The result is that the 
wave function, eq.(35), does not have a well defined limit as t — > 0~. In terms of expectation 
values, this divergence result in the expectation value for P blowing up. One finds that, 

< P >~ (-t) 2{1 - p) -+ oo. (39) 

We have considered the ground state wave function above. In the subsection that follows, 
we will give a general argument for why the same diverging phase factor arises for the wave 
function of any state. This general argument will also include quartic terms. 

3.3 General Analysis of Wave Function Near t — 

The behavior of wave function can be analyzed quite generally in the vicinity of t = 0. It is easy 
enough to carry out the analysis in the the full quantum mechanics system, eq.(7), including 
the quartic interaction terms. 

The Schrodinger equation takes the form, 

~Y d ^ + e ~* v (^ = ld ^> ( 40 ) 

where the potential, V(x), is 

V(x) = ^lu 2 x 2 + x 4 . (41) 



9 For p = 1, the divergence goes like log(— t). This follows from standard properties of Bessel functions, and 
also from the general discussion in the next subsection. 
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Since e* vanishes near the singularity let us begin by assuming that the potential energy term 
on the lhs of eq.(40) dominates, we will verify below that this assumption is self-consistently 
true. This gives, 

e-*V(x)ip ~ idti/j, (42) 

which can be easily solved to give, 

^(x,£) = e- iG(t)y(£ Vo(;r), (43) 

where, 

G(t)= jdte-* = - { j^, (44) 

and ipo(x) is a time independent integration constant. Here we have used eq.(ll) for the dilaton. 
For p > 1, we see that G(t) diverges at the singularity, leading to a diverging phase factor in the 
wave function. This divergence is a general feature, independent of the initial state, ipo(x).In 
the quadratic case where 

V(x) = \^i\ (45) 

we see that this phase factor agrees with what was obtained in the exact solution for the ground 
state, eq.(38). 

To check the self consistency of our assumptions let us evaluate the contribution due to the 
kinetic energy term in eq.(40) on the solution, eq.(43). It is useful to analyze the two cases 
p > 1, and p < 1 separately 10 , in both cases we see below that the kinetic energy term is 
subdominant compared to the potential energy term. 

For p > 1 the leading contribution to the kinetic energy comes when the spatial derivatives 
act on the phase factor, and not on ipo(x). This gives, 

- y^ * y [G(t) 2 V\xf + iG(t)V"ty. (46) 

Since G(t) diverges the dominant contribution comes from the first term on the rhs leading to 

-yfliV' ~ y [G(t) 2 V'(x)^ ~ t 2 e-*V'(xf^ (47) 

where we have used the behavior for G(t) in eq.(44). Comparison with the potential energy 
term in eq.(40) shows that this contribution is suppressed by an extra power of t 2 . 

For p < 1 the leading contribution comes when the spatial derivatives act on ipo{ x )- This 
gives, 

- yflgV' - -ye-* G ^Vo(z)". (48) 



10 A similar analysis can also be carried out when p = 1. 
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Comparing with the potential energy term in eq.(40) we see that this term is suppressed by an 
extra power of t 2p . 

From the point of view of the bulk dual cosmology, we are especially interested in the 
question of whether the state can be continued past t — 0, with the dilaton varying for example 
like, 

e* = |t| p . (49) 

We have found above that the wave function for a general state in the quantum mechanics model 
does not have a well defined limit as t — »■ CT, when p > 1. This means it is not meaningful to 
ask about its continuation for t > in this case. To obtain this continuation one would need 
to impose that the wave function at t — is continuous, i.e., meets the condition, 

ip(x, t = CT) = ip{x, t = + ). (50) 

This condition cannot be imposed if Lim t _ >0 -ip(x,t) does not exist. 

3.4 The Energy Blows up at t — > 

We continue with our general analysis of the wave function in the quantum mechanics system 
in this subsection and find that for a generic state, and all values of p > 0, the energy at t — »■ CT 
diverges. We will work below with a general potential V(x). 

The Hamiltonian operator, H, is given by the left hand side of the Schrodinger equation, 
eq.(40), eq.(59). We have argued above that the kinetic energy contribution is subdominant to 
the potential energy, near t — 0, so that, 

< H >~ e~* < V > . (51) 

The expectation value of the potential, 

< V >= [ dxV(x)ip*(x,t)ip(x,t). (52) 



Substituting for the wave function from eq.(43) we see that the phase factor drops out so that 
< V > near t — is given by 



<V>= / dxV(x)\i/jo(x)\, (53) 

and is time independent. This means the leading time dependence in < H > comes from the 
prefactor e~* in front in eq.(51), leading to the conclusion that the energy diverges as 

<H>^ (-t)~ p (54) 

when t — > 0. Note that this conclusion holds for all p > 0. 
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This conclusion can be avoided if the state is such that < V > vanishes. This issue is 
examined further in some detail in Appendix C. The conclusion, after analyzing sub leading 
corrections which could also have been potentially divergent contributions, is the following: 
Unless p > 2, in which case a divergent contribution to the energy arises from the kinetic 
energy term, the subdominant contributions to the energy do not diverge as t — ► 0. Thus, 
requiring that < V > vanishes is enough to ensure that the energy stays finite. 

Now even if < V > vanishes we should note that the expectation value of < H 2 > will 
diverge. From the discussion above it follows that 

< H 2 >~ e^ < V 2 > . (55) 

In general in a state where < V > vanishes, < V 2 > will not vanish. This means that even 
in those special states where the expectation value of the energy stays finite as t — ► 0, the 
fluctuations about this finite mean value will diverge n 

4 Analysis in Field Theory 

In the previous section we have analyzed the quantum mechanical model, eq.(7) extensively. 
Here we return to field theory, first discussing the toy model field theory, eq.(6), and then 
turning to the deformed M = 4 field theory eq.(3). 

The lessons from the study of the quantum mechanics model can be directly applied to the 
field theory, eq.(6). Carrying out the field redefinition in the full field theory gives rise to a 
Lagrangian, upto a surface term, 

L = --{dX) 2 -m 2 {t)X 2 -e*X\ (56) 

with m 2 being a tachyonic time dependent mass, 



n 



2 



m 2 = ~ (57) 



with, 

a 2 = 

4 



a 2 = % + 2), (58) 



which diverges as t — ► 0. 

The first lesson which carries over from quantum mechanics to field theory is that contrary 
to what one might have guessed at first, the quartic term cannot be neglected near t = 0. 



11 In the next section we will apply the discussion of this and the previous subsection to the M = 4 gauge 
theory. In that case, both < H >,< H 2 > should scale like N 2 - the number of colors. This means that the 
fluctuations in energy will be suppressed in the large N limit. This suggests that to leading order in 1/-/V, and 
for p < 2, the vanishing of < V > is sufficient to ensure that the expectation value of energy stays finite when 
t-»0. 
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The second lesson is that the variable X is not so convenient to work with and the analysis is 
more conveniently carried out in terms of the original variables X. The third lesson is that the 
analysis is conveniently carried out in terms of the Schrodinger picture. This last lesson is not 
easy to apply in field theory, since typically the Schrodinger picture has not been used in this 
context. Nevertheless with the experience of the quantum mechanics model in mind we will 
in this subsection analyse the field theory in the Schrodinger picture; various caveats will be 
discussed in the next subsection. 

In the Schrodinger picture in field theory the state of the system is described by a time 
dependent wave functional, ip[X(x),t], which satisfies the Schrodinger equation, 

-\jd 3 x e"^t + e-*V[Xty = W- (59) 

Here the potential energy, V[-X"], is a functional given by, 

V[X] = Jd 3 x{^(d l X) 2 + X 4 }. (60) 

The first term on the l.h.s. in eq.(59) is the kinetic energy, and the second term is the 
potential energy. We see that, like in quantum mechanics, the kinetic energy term has a 
prefactor e*, while the potential energy has the prefactor, e~*. This suggests that the kinetic 
energy term is once again subdominant close to the singularity, leading to the solution, 

il)[X{J), t] = e-^ G ^ v ™^ [X(x% (61) 

We see that the wave functional has a phase factor which diverges, as in the quantum mechanics 
case, resulting in a singular limit for the wave function if p > 1. Moreover with the kinetic 
energy being subdominant, near t — 0, the Hamiltonian is well approximanted by, 

< H >~ e _# < V > . (62) 

As is the quantum mechanics model the phase factor drops out in the expectation value of 
< V > near t — leading to the conclusion that in field theory as well, < H > goes like, 
eq.(54), and therefore blows up for all p > 0. 

The central assumption here is that the kinetic energy is subdominant compared to the 
potential energy near t = 0. This was shown to be self consistently true in the case of quantum 
mechanics. We will analyse this issue for the field theory in the next subsection. 

Before proceeding let us also mention that we analyze the quadratic theory in further detail 
in appendix D. We consider a dilaton profile of the form, e* = \t\ p , and evolve the field theory 
in this background, starting from the vacuum into the far future. It is useful for this purpose 
to regulate the dilaton profile near t = in a manner we make more precise in the appendix. 
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Our conclusion is that particle production always occurs. For p < 1 the particle production is 
finite, while for p > 1 it becomes infinite as the regulator is taken to zero. 

Turning now to the gauge theory, eq.(3), we see that the coupling of the dilaton to the 
quartic scalar potential and the fermionic Yukawa terms are proportional to positive powers 
of e* and can be neglected when the dilaton is very small. The Fermions and scalars have 
canonical kinetic terms. The gauge field in contrast has a non-canonical kinetic energy term, it 
is the analogue of the X field in the toy model, eq.(6). As t — ► 0, and the dilaton vanishes, it is 
this gauge kinetic energy term which will determine the behavior of the system. Accordingly, 
in the analysis below we focus on the pure gauge theory, without fermions and scalars, and 
with action, 

S = f^x(~)TrF ltv F^. (63) 



The equation of motion is, 
We work in Coulomb gauge, where, 



D lt {e- 9 F< a ') = 0. (64) 



A = 0. (65) 

Consider now the non-interacting theory. The equation of motion, eq.(64), for, v = 0, 
becomes, (this is the Gauss Law constraint), 

d Q {&Ai) = 0. (66) 

Thus the longitudinal part of the gauge field is time independent. We can now do an additional 
time independent gauge transformation to set 

djA j = 0. (67) 

The equation of motion, eq.(64), with v — i, then become, 

9o(e-*9 ^) + e^dj&A* = 0, (68) 

for the two transverse components satisfying, eq.(67). Eq.(68) is exactly the equation we have 
for a scalar field with Lagrangian, 

L = e-*(dX) 2 . (69) 

Thus at the quadratic level, the analysis for the gauge field reduces to that of the scalar field 
considered above. There are two transverse components coming from each gauge field. We have 
neglected the color degrees of freedom above. They are easily incorporated and give (N 2 — 1) 
degrees of freedom for each of the two transverse components. 

Interactions give rise to cubic and quartic terms. In Coulomb gauge, these terms do not 
depend on any time derivatives. Thus, in the Hamiltonian they only contribute to the Potential 
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energy and not to the Kinetic energy. The interactions can therefore be included in a way very 
similar to the quartic terms in the toy model for X. 

In particular, we are interested in the wave function in the Schrodinger picture. In this 
picture the operators are time independent. The total potential energy is given by energy due 
to the magnetic field, 

V[A(x)] = J <Px±Tr(F ij F i t). (70) 

Motivated by the quantum mechanics model on the previous section, and as in the scalar 
field theory above, we now take the potential energy to dominate in the Schrodinger equation 
near the singularity. As a result the wave function has a phase given by, 

if>[Ai(x),t] = e-^W^MMx)}. (71) 

The phase factor is identical to that found in eq.(43) and diverges, as t — ► 0, if p > 1, resulting 
in the wave function being singular at t — > 0. The wave function above can be regarded as 
being dependent on only the transverse components of the vector potential. Alternatively, we 
can take the wave function to be dependent on a general gauge potential (with A — 0), and 
then impose Gauss' law, 

*& = "■ < 72 > 

on it. 

Similarly one can calculate the expectation value of the energy near t — 0. It has the same 
form as in eq.(54), with the expectation value of the potential energy now being given by, 

< V >= J DAilMM^WViAiix)], (73) 

where V[A^(a;)] is given in eq.(70). We see that the energy diverges in the gauge theory as well, 
as (— t)(~ p \ for all p, unless the state is such that < V > vanishes. 

We conclude this section with two important comments. First, naively one would have 
thought that as the dilaton becomes small perturbation theory should become a good approx- 
imation. However we have seen in our analysis of the toy model in the last section that this 
is in fact not true. In the X description for the toy model, the time dependent dilaton drives 
the system to large values of X resulting in the quartic term being non-negligible. A similar 
argument also holds in the gauge theory. The cubic and quartic interactions terms are not small 
near t — 0, and as a result perturbation theory is not a good approximation. Secondly, it should 
be emphasized that our analysis is valid for finite N and finite g\ M X . In fact from eq.(71) we 
see that the behavior of the wavefunctional near t — is essentially determined by G(t) and 
is independent of the details of the potential V[Aj(a;)]. In the dual closed string theory, this 
means that this conclusion is valid in the presence of string and quantum corrections. 
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4.1 A more critical look at the Field Theory Analysis 

The central assumption in the field theory discussion above was that in the Schrodinger equation 
the potential energy which scales like e~* dominates over the kinetic energy, which has a 
prefactor e* in front of it. This assumption was motivated by our earlier analysis in quantum 
mechanics. However, field theory differs from quantum mechanics in having an infinite number 
of degrees of freedom, and one might worry that this introduces additional subtleties and 
complications 12 . We turn to an examination of these issues below. 

The wave function for the ground state of the Harmonic oscillator with action eq.(14) in the 
x description is given by eq.(35). This is an exact result. If the potential energy is dominant 
the wave function has the form, eq.(43). From the exact result we can ask how close to t — 
must one come for this approximation to become a good one. The phase factor in the wave 
function eq.(35) contains the factor, (4)* + *. As discussed in eq.(150) of Appendix B this 
phase factor has a power series expansion in tuo, near t — 0, of form, 

The ellipses stand for terms which are supressed when 



2c 2 uj 2 t + ■■■. (74) 



tu> < 1. (75) 

Now the first term in eq.(74) vanishes due to eq.(17). The second term is therefore the leading 
one and it is easy to see that this gives agreement with eq.(43). 

The conclusion is that the potential energy dominates over the kinetic energy term when t 
is small enough and satisfies the condition in eq.(75). We can also understand this on the basis 
of the general arguments given in section 3 for the self consistency of this approximation. From 
eq.(47) we see that the KE term goes like, 

-yfl|V ~ tVM^V ~ i W 2 S)V, (76) 

while the potential energy term is, 

e-*V(x)ij; ~ e- $ (cu £)V. (77) 

Thus for the latter to dominate, eq.(75) must be true. 

Consider now a field theory in the quadratic approximation, since there are an infinite 
number of modes, for any non-zero and arbitrarily small time t, there will always be some 
modes with high enough frequency for which the condition eq.(75) is not met and thus for 
which the wave functional will not be well approximated by eq.(61), or eq.(71). Including 



12 We are grateful to David Gross in particular for emphasising this point to us. 
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interactions will couple these modes to low-momentum frequency potentially making the wave 
functional for all modes to be different from eq.(61), eq.(71). 

It is of course well known that in dealing with the infinite numbers of degrees of freedom 
in a field theory it is first useful to introduce a cut-off or regulator, which makes the number 
of degrees of freedom finite and then ask what happens as the cut-off is removed. For ease 
of discussion consider a momentum space cut-off A ( in the gauge theory one needs a more 
sophisticated regulator to preserve gauge invariance, but this will not change the essential 
points in our discussion). The above analysis suggests that if we take the time t to vanish while 
keeping A fixed and finite, then the potential energy should dominate for times t meeting the 
condition 

£A < 1 (78) 

and our conclusions in the previous section will be correct near t — 0. More generally one might 
expect that these conclusions are valid as long as we take t — > 0, before we take, A — »■ oo. 

Unfortunately, it is not easy to make this argument precise. Additional complications can 
arise in field theory due to the effects of renormalisation 13 . Usually renormalisation in field 
theory is discussed in terms of an effective Lagrangian which changes under RG flow. For our 
purposes in the discussion above the Schrodinger picture has been more useful. Including the 
effects of renormalisation in the Schrodinger picture though is a complicated issue that we have 
not fully sorted out. Presumably the Hamiltonian which governs the time evolution of the wave 
functional needs to be made well defined by appropriate operator ordering and this introduces 
the effects of renormalisation. 

One would expect that at least some of the consequences of renormalisation can be incor- 
porated by first constructing an effective Lagrangian by integrating out high frequency modes 
and then using this effective Lagrangian to construct a Hamiltonian that governs the time evo- 
lution of the surviving low-frequency modes. Since the M = 4 theory is conformally invariant 
any renormalisation of the effective Lagrangian must be due to the time dependence of the 
dilaton and thus operators which are induced by this renormalisation must have coefficients 
proportional to derivatives of the dilaton. In turn such operators could then also change the 
Hamiltonian resulting in extra operators in it with coefficients proportional to time derivatives 
of the dilaton. 

The essential reason why in our discussion the potential energy dominates is that it scales 
like 

V ~ e"*, 

whereas the kinetic energy term has a prefactor e* in front of it that supresses it. However 
suppose as an example of the consequences of renormalisation the potential energy acquires an 



13 We thank the referee for emphasising this point to us. 
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extra term which arises at one loop so that it now has the form, 

Here 0\ is the operator corresponding to the magnetic field energy, A is the cutoff scale and 
O2 is the additional operator which arises at one loop. If the condition 

p - 2 < (79) 

is met, this second term could get important close to t — 0. This will mean that one has to 
include additional loop effects that arise beyond one loop as well. Resumming these effects 
could well lead to a much smaller potential energy. For example, if these corrections take the 
form of a geometric series, we would get schematically 

we would get after resumming 

V~e~*[ . 1 ~ -^ , (80) 

L A 2 -e*($) 2J ($) 2 V ; 

sufficiently close to t — 0, if the condition eq.(79) is met. We see that the effects of renormali- 

sation can therefore supress the potential energy term. In particular if 

A 2 
e -^- ~ e 81 

$2 v ; 

this supression would make the the Kinetic energy term comparable. It is easy to see that 

eq.(81) will be met for small enough time,t, if, 

p < 2/3. (82) 

The summary is that our analysis in the field theory is not complete and our conclusions 
about the gauge theory being singular should be taken as being suggestive but not conclusive. 
To analyse the gauge theory in a well-defined manner one must introduce a regulator. Once 
this is done a very rich set of counter terms are allowed in the process of renormalisation, and 
such counterterms can potentially invalidate our conclusions. This would happen if they make 
the potential energy comparable to the kinetic energy or even smaller than it near t — thereby 
significantly changing the form of the wave function and potentially making the gauge theory 
non-singular. Whether this happens or not requires a detailed understanding of renormalisation 
in the gauge theory in these time dependent background. This is a fairly complicated subject 
and we leave it for the future. 

We end with some comments. Introducing a UV regulator in the gauge theory is dual to 
introducing a boundary in the bulk that regulates the IR behaviour. Some consequences of 
renormalisation, in the null dependent case, have been worked out in [5], where indeed terms 
in the Wilsonian action with coefficients proportional to i -^- were found. 
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5 Cosmological Solutions 

The motivation for this investigation came from trying to understand some cosmological so- 
lutions [1, 2, 3]. These solutions can be thought of as deformations of AdS§ x S 5 and have a 
dual description in terms of the M = 4 theory with a time dependent dilaton. As the dilaton 
becomes small on the boundary the bulk curvature becomes larger and larger, eventually be- 
coming singular at t — 0. In this section we ask what the above analysis in the gauge theory 
teaches us about these cosmological solutions. 

We begin with a brief review of these solutions and then return to the gauge theory later. 

5.1 The Gravity Solutions 

The solutions arise in IIB theory and are deformations of AdS$ x S 5 . The S 5 factor is unchanged 
in the deformations, and accordingly we will omit it below and only discuss the solution in the 
remaining five dimensions. 

The first solution we consider has the 5-dimensional metric, 

1 dr 2 

ds 2 = A,[dz 2 + \Sinh{2t)\\-dt 2 + -— - + r 2 [dQ 2 + sin 2 0# 2 )]], (83) 

z 2 1 + r 2 

and dilaton, 

e^ = g s \tanht\^. (84) 

This solution was discussed in [3]. 

In the far past, as t — ► — oo, the dilaton goes to a constant, and the metric becomes AdS§. 
One can see this by going to coordinates, 

R 



, ,_ e"' = vV - R 2 , (85) 

in which the metric, eq.(83), and dilaton, eq.(84), take the form, 

ds 2 = ^[dz 2 + |1 - 2 ^ R2)2 \[-dr] 2 + dR 2 + R 2 dQ 2 2 ]]. (86) 

The far past, t — > — oo, corresponds to (rf — R 2 ) — ► oo, it is clear from eq.(86) that the metric 
asymptotes to AdS$ in this limit. The dilaton in these coordinates is, 

y-fl 2 -i |V3 (87) 

\rj2_ R 2 + 1 \ ■ y*<) 

At t = the solution, eq.(83) has a singularity. The curvature scalar diverges like R ~ A- as 

The dilaton vanishes as t — > 0, eq.(84). Thus the singularity occurs at weak string coupling. 
This singularity is the main focus of our analysis. 
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The region t < 0, which is the region of spacetime before the singularity, maps to n 2 — R 2 > 1, 
in the (77, i?) coordinates, while the singularity, which is at t — 0, maps to the locus, rf — R 2 = 1. 
Another 5-dim solution is given by 

rfs 2 = \[dz 2 + \2t\[-dt 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 }, (88) 
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with dilaton, 



e *(t) = ^| f |V3. (89) 



This solution does not asymptote to ArfS'5 in the far past, as t — > —00. However its behavior 
at the singularity, as t — > 0, is very similar to the solution discussed above, eq.(83), eq.(84). 
The z = const hypersurfaces in both metrics are of the FRW form. The difference between 
the two metrics is that this 4 dim. FRW cosmology has constant negative curvature in the 
first case, eq.(83) while it is flat in the second case, eq.(88). This difference is increasingly 
unimportant near the singularity, where the dominant source of stress energy is provided by 
the diverging time derivative of the dilaton, rather than the spatial curvature. Since the dilaton 
is essentially identical near the singularity, at t — ► 0, in both cases, the resulting spacetimes 
also are essentially the same. 

In section 6 we will explore some additional cosmological solutions and comment on their 
gauge theory duals. Some of these solutions differ from the two solutions discussed above at 
early times but their behavior near t — becomes the same as in the solutions above. 

5.2 The Gauge Theory Duals 

For purposes of studying the Field theory dual, we start with the first bulk solution considered 
above, eq.(83), eq.(84), or equivalently eq.(86), eq.(87). This solution asymptotes in the far 
past to AdS§ with a constant dilaton. As discussed in [1], [2], this corresponds to starting in 
the far past with the N = 4 SYM theory in the vacuum state. The spacetime, eq.(83), has a 
boundary at z — ► 0. We see that the metric on the boundary is conformal to flat space. As was 
discussed in [3] , we can then take the metric of the spacetime in which the dual gauge theory 
lives to be the 4 dim. Minkowski-space metric. 

Since this is an important point let us pause to briefly comment on it further. In general a 
Weyl transformation in the boundary theory corresponds to a Penrose-Brown-Henneaux (PBH) 
transformation in the bulk. The explicit PBH transformation which gives rise to a flat boundary 
metric for eq.(83) was found in [3]. The metric in eq.(83) has a second order pole as z — ► 0. 
The metric on the boundary defined by 

ds\ = Lim z ^ z 2 g )MU dx tl dx u , (90) 
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where x^ denotes coordinates on the boundary, is, 

civ 

ds\ = \Sinh{t)\[-dt 2 + y- - + r 2 {dd 2 + sin 2 6d<f) 2 )]. (91) 

After the PBH transformation the resulting 4 dimensional metric is given by, 



ds\ = e' 2t [-dt 2 + -^— + r 2 (d9 2 + sin 2 9d<f) 2 )}. (92) 

1 + r 

This is in fact flat space in Milne coordinates. It is easy to see this. The coordinate transfor- 
mation, eq.(85), turns this metric into the familiar Minkowski metric, 

ds\ = -drj 2 + dR 2 + R 2 dQ 2 (93) 

The region to the past of the singularity, is given by, — oo < t < 0, in eq. (92), and maps to 
the region rj 2 — R 2 > 1, which is part of one of the Milne wedges. The rest of this Milne wedge 
is given by < t < oo, to which the metric eq.(92) automatically extends. The boundary of 
the Milne wedge lies at t — > oo. Starting from t — ► — oo one arrives at the singularity, at t — 0, 
before reaching the boundary of the Milne wedge. 

The dual gauge theory knows about the time dependence of the bulk through the varying 
dilaton. The exponential of the dilaton is equal to the coupling constant of the Yang Mills 
theory, eq.(2). Thus the varying dilaton gives rise to a varying Yang Mills coupling. The 
dilaton depends on the Milne time coordinate, t, eq.(92), and takes the form, eq.(84). 

In summary we see that the dual gauge theory to this cosmological solution lives in flat 
space with a varying dilaton which vanishes at the singularity. The analysis in the preceding 
sections can be now be used to determine the nature of this singularity. We turn to this in the 
next subsection below. 

Before that, it is also useful to discuss the dual to the second cosmological solution intro- 
duced in the previous subsection, with metric and dilaton given by, eq.(88), eq.(89). In this 
case the dual gauge theory also lives in flat space, but the dilaton depends on Minkowski time 
instead of Milne time. To see this note that from eq.(88) it follows that the boundary metric 
as defined by eq.(90) is conformally flat. After a suitable PBH transformation the boundary 
metric becomes that of flat Minkowski space, 

ds 2 = -dt 2 + (dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 - (94) 

The dilaton dependence is given in eq.(89), we see that it is only a function of Minkowski time. 
The first solution, eq.(84), has a dilaton which goes to a constant in the far past, and the 
dual gauge theory starts in the vacuum state of the N = 4 theory as t — »■ — oo. In contrast 
in the second solution, eq.(89), the dilaton blows up in the far past, this makes the dual map 
to the boundary theory less clear. It will turn out that the behavior at the singularity of the 
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two theories is similar 14 , near the singularity, and a little easier to analyze in the second case, 
where the variation is with respect to Minkowski time. For this reason, and for the limited 
purpose of asking questions near the singularity, we will focus some of the following discussion 
on the second hologram. 

5.3 Gauge Theory and Gravity 

In this section we relate what was learnt in the analysis of the gauge theory above to the specific 
cosmological solutions of interest. 

We see from eq.(84), eq.(89), that the solutions correspond to 



P 



V3 (95) 



In particular this means that for these solutions, p > 1. We now see that, with the provisos 
discussed above, the analysis in the gauge theory suggests that the system becomes genuinely 
sick as t — ► 0~, . In particular the wave function in the Schrodinger picture, acquires a wildly 
oscillating phase, eq.(43), and thus does not have a good limit, as t — ► 0. This also means the 
state cannot be sensibly continued past t — 0. We remind the reader that these conclusions are 
not definitive, in particular various caveats discussed in section 4.1 apply here, and there are 
additional issues, having to do with renormalisation that we have not adequately discussed in 
this work. 

If true, the conclusions should hold regardless of the state of the system. In the far past, 
for the solution, eq.(83), eq.(84), the state is known to be the vacuum of the N = 4 theory. 
However, as time progresses and the dilaton becomes smaller the state evolves. The analysis in 
the gauge theory leading to the wave function of form, eq.(71), is only valid very close tot — 0, 
where the state would be different in general from the vacuum. 

It also follows from the gauge theory, eq.(54), that the energy diverges like, (— £) 3 , as 
t — > 0. If our approximations hold, the only way to avoid this conclusion would be if the 
system, which starts in the vacuum state in the far past, evolves to a non-generic state near 
t = for which < V >, eq.(52), vanishes. While this seems unlikely, since such states are non- 
generic, one cannot rule out this possibility. Note however that even in this case the discussion 
of the previous two paragraphs would continue to hold and the system would become sick as 
t — *■ 0~. Also, as was mentioned in section 3, even in such a non-generic state the fluctuations 
in the energy, which are suppressed at large iV would still diverge and the gauge theory at finite 
N would still be singular. 

The cosmological solution in eq.(83), eq.(84) maps to a boundary theory where the dilaton 
depends on Milne time, rather than Minkowski time. This does not make an essential difference 



14 This similarity is in parallel with the fact, mentioned in the previous section, that the two corresponding 
bulk solutions also behave similarly near the singularity. 
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to the analysis in the gauge theory which in the above sections was carried out for the dilaton 
being a function of Minkowski time. This issue is analyzed further in the appendix E. The 
Milne case correspond to having a non-trivial metric in the boundary theory. The dilaton 
vanishes at the point t — which is perfectly smooth in Milne coordinates, thus the nontrivial 
metric does not make an essential difference to the discussion of the singularity. 

More generally one can consider other cosmological solutions, which are different in the far 
past, but which also behave like the two examples discussed above, near t — where the dilaton 
vanishes. In all these cases, as long as the metric on the boundary is well behaved at t — 
the above discussion should apply. In section 6, we give examples of BKL cosmologies, with 
the dilaton, which at late times asymptote to the example, eq.(88), eq.(89). It follows from 
the discussion above that, subject to the caveats mentioned above, the singularity in all these 
solutions then is a genuine sickness of the theory. 

5.4 Concluding Comments 

We conclude this section with some more comments on the relation between the gravity solutions 
and their gauge theory description. 

For the supergravity solutions, eq.(83), eq.(84), eq.(88, eq.(89), the stress energy tensor was 
calculated in Section 5 of [3]. For small t the stress tensor diverges like, 

N 2 
T,u ~ ir - (96) 

This calculation was made in the super gravity approximation which breaks down when the 
dilaton becomes small enough. At very small values of the dilaton the gauge theory analysis 
carried out in section 4 becomes valid, subject to the various caveats discussed in section 4.1. 
We have seen in section 4 that the energy density according to this analysis goes like 15 

< p >~ iV 2 /|£| p . (97) 

(The factor of iV 2 was not explicitly displayed in the discussion in section 4. It is easy to 
see that in the quadratic approximation, where the gauge theory is free, the N 2 independent 
color degrees of freedom, in the presence of the time dependent dilaton, give rise to an energy 
density that scales in this way with N. In the presence of interactions, the energy density 
should continue to scale like N 2 , to leading order in N.). Since p = \^3 < 4 we see that the 
growth of energy in the regime where the dilaton is very small, is much slower than that in 
supergravity. 



15 The initial state in the far past is the vacuum which is translationally invariant. Also, the dilaton only 
depends on time and docs not break this translational symmetry. Thus the potential energy in cq.(70) scales 
like the volume, leading in turn from eq.(51) to an energy density which is finite and given by eq.(97. 
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In the cosmological solutions we are studying here, the energy being pumped in by the 
dilaton does not lead to the formation of a black hole in the regime where supergravity is 
valid, for t < 0. This means that in the dual gauge theory the energy being pumped in is not 
thermalizing. To understand this we note that in the strongly coupled gauge theory, which is 
dual to supergravity, the thermalization time scale r, is expected to be of order 

1 /V 2 

T ~-~(—)W (98) 

T P 

where p is the energy density, and T is the temperature. From stress energy tensor, eq.(96), 
we see that 16 

r ~ t. (99) 

Thermalization would occur if the external source is pumping in energy on a time scale much 
slower than r. In the situation at hand, this time scale is of order, 

p 

tpump ~ -■ (100) 

Using eq.(96) we see that t pump is also of order t. Thus the rate at which the dilaton is 
pumping in energy equals the thermalization or relaxation rate of the system. This explains 
why thermalization does not happen and in the dual a black hole does not form. Another way 
to see this is that from eq.(98), eq.(99) it follows that, 

^~0(1), (101) 

so that the temperature changes too rapidly for thermalization. 

We have seen that within our approximations, a dilaton which decreases all the way to 
zero results in a singular gauge theory. Effects of renormalization could, in principle, tame the 
singularity. In any case, if the dilaton profile is modified such that e* never vanishes, but can 
beome small at t — 0, a smooth time evolution beyond this time is possible. Let us consider 
a situation where the dilaton varies in a smooth fashion reaching a minimum value and then 
increasing again, approaching a constant in the far future which corresponds to a large value for 
the 'tHooft coupling. The time dependent dilaton will typically lead to the system having some 
non-zero energy density in the far future, and one would expect that given enough time this 
energy would thermalize. This suggests that in the dual closed string description supergravity 
will eventually become a good approximation in the far future, and the dual geometry in the 
far future will be that of a black hole in AdS^. 

An interesting question to ask is how the formation of the black hole depends on the dilaton 
time variation. In particular whether the black hole is always a good description of the geometry 



16 Note that eq.(96) is valid when \t\ <C 1 for the solution eq.(83). This is consistent with the supergravity 
approximation being valid, as follows from eq.(104). 
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by the time the supergravity approximation becomes good in the future, or whether for a 
suitable dilaton time profile the formation of the black hole can be delayed, till much after the 
supergravity description becomes valid. This question is difficult to answer with our current 
level of understanding. In particular, in the course of the time evolution of the dilaton it goes 
through a region where the 'tHooft coupling is of order unity. This region is very difficult to 
analyze, since neither the gravity nor the gauge theory descriptions are tractable then. 

Finally, we have not been very precise in the discussion above about exactly when the 
supergravity approximation breaks down. The metric for the solutions we are considering is in 
eq.(83), eq.(88). For small t the curvature goes like, 

z 2 

(-t) 3 

In this expression we have set 

RA d s = {^g s N) l/A (103) 

to unity, and we are also working in 5 dim. Einstein frame. After taking this into account one 
finds that the string frame curvature is given by, 



strin l Z<2 



U3-^) m 



As t — ► 0~ we see that the curvature blows up, resulting in a singularity as we have discussed 
above. The supergravity approximation breaks down when the curvature gets to be of order the 
string scale. We see that when this condition becomes valid depends on the radial coordinate 
z. In particular, for all finite t, z — > oo has diverging curvature. This point corresponds to the 
locus where the past and future horizons of AdS§ (in the Poincare coordinates we are using) 
intersect. The significance of this curvature singularity is unclear to us 17 . As time increases, 
the region with high curvature (of order string scale or more) grows, moving to smaller z, 
eventually leading to a singularity for all z as t — ► 0. 

The calculation of curvature, eq.(104) makes it clear that the singularity (in string frame) 
is due to two effects which are tied together in these cosmological solutions. First, the dilaton 
vanishes causing the string frame curvature to blow up. Second, an infinite amount of energy 
is dumped into the system due to the vanishing dilaton source. If the dilaton profile is altered 
so that it attains a non- vanishing minimum value at t = 0, the total energy put into the 
system would be finite. If the minimum value of the dilaton is small enough one expects that 
the curvature in string units at t — still becomes much bigger than unity, resulting in the 
breakdown of supergravity. 



17 Note that the past Poincare horizon at z — * oo, t — * — oo, with \z/t\ held fixed is non-singular, as was 
discussed in [3]. 
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6 Behavior Near Singularities and BKL Cosmologies 

We saw in section 5 that the two solutions which were considered had very similar behavior near 
the t — singularity. This was because the stress energy near the singularity was dominated 
by the dilaton which had essentially the same behavior near t = in these two cases. 

The behavior of the dilaton near the singularity is in fact shared by a much larger class of 
cosmological solutions. The key point is that the Einstein frame metric for the class of solutions 
we have considered may be written in suitable coordinates as [1] 

1 



ds 2 



z 2 



dz 2 + g ia ,{x)d3?dx v +dVtl (105) 



This solves the 10 dimensional supergravity equations provided 

Rp, - ^V M $V,$ = 0, 



Kd^y/HJ gVd v $) = 0. (106) 



V-9 

and the five form field strength is standard 

F = uj 5 + *w 5 (107) 

In other words, any solution of 3+1 dimensional dilaton gravity may be lifted to a solution of 
10 dimensional supergravity. This means that we can use the well known analysis of Belinski, 
Lifshitz and Khlalatnikov (BKL) and subsequent work [13]- [18] to make useful statements about 
AdS cosmologies. 

For instance, we can consider a cosmological solution where the spatial 3-metric is one of 
the general homogenous spaces in the Bianchi classification (see e.g. [13] for a lucid treatment 
in the 4-dim context), with vanishing go yCe components: 

ds\ = -dt 2 + r] ab (t)(e a a dx a )(e b pdxP) , (108) 

where (e^,dx a ) are a triad of 1-forms defining symmetry directions 18 . i] a b(t) are general time- 
dependent coefficients which can be solved for from the Einstein equations, with components 
decomposed along the frame. Assuming a spatially homogenous dilaton gives d a Q = e"9 a $ = 0, 
with <9o$ nonvanishing, so that R? a -, vanish, with Rq = |(<9o$) 2 . More details on the Bianchi 
IX solution can be found in Appendix F. 

The main point of BKL is that close to a spacelike singularity, physics becomes ultralocal. 
For dilaton driven cosmologies, this results in a Kasner-like solution in which the time dependent 
part of the dilaton is precisely of the form we have been analyzing. 



18 Starting with the 1-form triad (e^dx a ), a labelling the vectors in the triad, we can obtain the dual vectors 
e", satisfying e a .ei> — 5%. Then the symmetry algebra acting on the homogenous space (i.e. the spatial metric) 
in question is obtained as the algebra of the differential operators X a — e"d a . 
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A simple illustration is provided by a general conformally flat boundary metric and dilaton 
ds 2 = F(x, t) i-dt 2 + dx 2 ) $ = $(x, t). (109) 



Near a singularity (which may be chosen to be at t = without loss of generality) we will assume 
that space derivatives may be ignored compared to time derivatives (which would typically blow 
up). However terms which contain mixed derivatives need to be retained [13]. This results in 
the following system of equations for $ and / = log F 

d t 2 f(x,t) + [d t f(x,t)] 2 = 

d t f(x,t)dMx,t) + d t 2 <S>(x,t) = 

3d t 2 f(x,t) + [dMx,t)} 2 = 

d t f(x,t)d l f(x,t)-2 d t dif(x,t)-dMx,t)dMx,t) = 0. (110) 

The general solution of (110)is given by 

$(x,t) = V3 ln(t + D(x)) + d(x), f(x,t)=\n(t + D(x)) + C 2 (x). (Ill) 

The last equation in (110) imposes the following relation on C 2 and C±, 

C 2 (x) = V3d(x). (112) 

By choosing D = 0, one can see that the behavior of the fields is the same as in eq.(88) and 
eq.(89) for the symmetric Kasner case. 

A similar result holds for a general class of diagonal metrics with non- vanishing Weyl tensor 
[15]. 

In fact, in dilaton driven cosmology for any space-time dimension greater than 3, the ap- 
proach to a singularity is characterized by a finite number of oscillations between Kasner-like 
solutions. Consider for example homogeneous cosmologies of type Bianchi IX 

ds 2 = —dt 2 + (a 2 (t)l a l/3 + alfymainp + a\{t)n a np) dx a dx t3 , (H3) 

where I, m, n are the three frame vectors e 1 , e 2 , e 3 . The Kasner-like solutions are obtained when 



5 ^ 5 



the spatial curvatures can be ignored, 

ai(t)~t p > ,$~alog(t) (114) 

where the Kasner exponents satisfy 

a 2 

Ep* = 1 > J2p" = 1 -y- (U5) 

i i 
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The effects of the spatial curvature results in oscillations between different sets of pis till all 
the pj's are positive. 

The transition between different Kasner regimes lead to an interesting attractor behavior. 
Consider the case of 3+1 dimensions Let p- denote a negative Kasner exponent and p + > 
being either of the other two positive exponents. These transitions can be expressed as the 
iterative map 

(n) (n) . (n) 

(n+1) -P- (n+1) P+ +2P- ®n n , „x 

"■ = T7^' p ' = 77^' """^ITW (116) 

for the bounce from the (ra)-th to the (n + l)-th Kasner regime with exponents Pi, pj. With each 
bounce, we see that a increases, shrinking the allowed space of {pi}, as detailed in Appendix 
F. For a > 1, the allowed window of {pi} pinches sufficiently forcing all pi > 0, at which point 
oscillations cease and the system settles in the pi > attractor region. The rate of increase 
of a is small for small a, since a n+ i — a n = a n ( ~ ~ 2 P ~ ), so that a system with near constant 
dilaton (a ~ 0) takes a long iteration time to "flow" towards the pi > attractor region. 

Furthermore many distinct initial Kasner regimes can flow to the same attractor point 
characterized by a set of positive pi's. The details of the derivation are given in Appendix F. 
Note that in the absence of a dilaton the number of such oscillations is infinite even though 
the proper time to the singularity is finite [13, 16, 15], simply because in this case all the pis 
cannot be positive. 

Finally, these attractor flows exhibit some degree of chaotic behavior, in the sense that small 
changes to the initial conditions give rise to drastic changes in the final endpoints, as elaborated 
in Appendix F. A quick glimpse at this is obtained by considering exponents {pi,P2,P3, a = 0} 
corresponding to a non-dilatonic asymmetric Kasner cosmology which oscillates indefinitely, and 
perturbing infinitesimally. Now a 2 = 2(1 — J2p'i ) is generically nonzero (although small). This 
latter set {p' { , a^O} thus flows to the attractor region, while the former oscillates indefinitely, 
showing that a small change in the former gives a drastically different endpoint. 

Note that in our setting we have frozen the 5 form field strength and all the other super- 
gravity fields. This is because we want to embed BKL type cosmologies in the AdS setup 
used in this paper. In general supergravity theories, a BKL type analysis shows that a general 
solution (which excites all fields) in the supergravities which follow from string theories and 
lid supergravity exhibit an infinite number of oscillations between different Kasner regimes 
[17, 18], similar to pure gravity [13, 16, 15]. 

The fact that the general BKL analysis for gravity-dilaton system can carried over to a 
discussion of a class of AdS cosmologies is interesting. However, the symmetric Kasner is the 
only solution whose Weyl curvature vanishes 19 . In the AdS context this means that it is only 



19 Analysing the Weyl tensor components shows that the Weyl tensor vanishes identically only for flat space 
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for this case that we can have PBH transformations to choose a flat boundary. In other cases, 
the metric on which the dual CFT lives is generically singular and one would expect that the 
gauge theory would be singular as well. 

A Dilaton Couplings in the Yang-Mills Theory 

In this appendix we discuss the Lagrangian which appears in eq.(3), especially the dependence 
of the dilaton in it. 

In the standard AdS/CFT dictionary the operator dual to the dilaton is determined by the 
superconformal symmetry to be an appropriate descendant of the chiral primary obtained by 
symmetrizing two scalars. Of particular importance to this paper is the fact that in eq.(3) the 
dilaton only couples to the kinetic energy term for the gauge fields and does not couple to the 
kinetic energy terms of the scalars and the fermions. To see this it is enough to consider the 
£7(1) theory. We follow the notation in [20]. The supercharges are Q^Qia- Here, the index / 
upstairs (downstairs) denotes a 4 (4) of SU(4), and a, a are indices for the two different spinor 
representations of 5*0(3, 1). 

The scalars transform like a real 6-dimensional representation of SU(4). For the limited 
purpose of carrying out the supersymmetry analysis it is useful to denote the scalar fields by, 
$[/ j], where the square brackets indicated antisymmetrization. The scalars satisfy the condition, 

($[„])• = \e IJKL ®KL. (117) 

The traceless symmetric product of two scalars gives rise to a field which transforms in the 
20 dimensional representation of SU(4). This is a chiral primary of the full superconformal 
algebra. We denote it by, 

Tukl = ®[ij]®[kl] ~ ^e/j^L$[PQ]$[R5]e PQi?5 - (118) 

Note this field satisfies the tracelessness condition, e IJKL TjjKL = 0. 

The operator which the complexified dilaton-axion couples to is, [20], 

= ^^QlQiQXQ^TjjKL. (119) 

The supersymmetry transformations in the U(l) theory are (upto possible numerical fac- 
tors): 

[QIQjk] = #*fc«-<&*Ja (120) 

{Q^Ia} = 5jF a p, (121) 



and the symmetric Kasner spacetime. For a generic asymmetric Kasner spacetime with exponents (pi,P2,P3), 
some of the nonvanishing Wcyl tensor components diverge as i^v-m-i ^ where p m is one of the exponents pi. 
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and, 

[#,^1 = 0. (122) 

Here \I/j a , F a ,$ denote the fermionic partners and the gauge fields. It is then easy to see 
that this gives 

6 ~ F 2 . (123) 

In particular O does not contain any coupling to the scalar or fermion kinetic terms. 

This result is consistent with eq.(3). It is also consistent with the statement that the dilaton 
couples to the on-shell Lagrangian once we allow for a total derivative term involving the scalars. 

A further check on eq.(3) may be obtained as follows. Consider deformed AdS$ x S 5 with 
a ten dimensional string frame metric g a p(x) and a dilaton $(x). The indices (a, f3) are 10 
dimensional indices, e.g. a = (fi, a) where // = 0, • • • 3 and a = 1 ■ • • 6. We will only consider 
backgrounds such that g ail = 0. Then the Lagrangian for the dual theory for small deformations 
is, 

C = v^ e"*Tr [-Ig^'^F^F^ 

~g^D,X a D v X b g ab + ±[X a ,X b }[X c , X d ]g ac9bd 

+ hr A e\ \-iD^ tt] + hr A e\ \X\ *] 9ab }. (124) 

Here g& denotes the determinant det(g IJ , 1/ ) and A denotes a frame index and e a A is the string 
frame vierbein. One way to see this is to consider the Yang-Mills theory in the Coulomb branch 
with SU(N) — ► SU(N — 1) x U(l). Then the effective action for the U(l) part should be given 
by the DBI action for a 3-brane in this geometry. In this action, the dilaton factor e* appears 
as an overall factor, provided everything is written in terms of the string frame metric. The 
leading order (two derivative) terms of this action can be obtained by simply replacing the 
SU(N) fields in the original YM action by the U(l) part. It is then easy to see that if we make 
this replacement in (124), we get the correct leading terms of the DBI action. 

We now need to express the Lagrangian (124) in terms of the ten dimensional Einstein frame 
metric G a p 

G a p = e-*' 2 g aP . (125) 

This leads to the Lagrangian 

C = y/=GTr[-±e-*G'*'G?" / F lu ,F lt , v , 

~G^D,X a D u X b G ab + ^e*[X a ,X b ][X c ,X d ]G ac G bd 

+I e »/«*r A {e E y A [-iDp, tf] + l -e^/^T A (e E ) A [X b , *}G ab ], (126) 
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where (e#)^ denotes the Einstein frame vierbein. Now consider the field redefinition for the 
fermion fields 

■q, ^ e -*/8tf. (127) 

This will absorb the dilaton factor in front of the quadratic term, but will give rise to an addi- 
tional term of the form 20 , ^Y A ^ [d^^i^E)^- However, ^T" 4 ^ = by virtue of the Majorana 
condition. 

In our setup, it is the Einstein metric which is flat. This gives rise to the equation (3). 

Finally, let us mention that eq.(3) is invariant under the conformal transformation, g^ v — > 
Q 2 g IIU , with the scalars and fermions transforming in the standard manner, once the \RX 2 
term is also included. 

B Reduction to Single Harmonic Oscillator and 
The Schro dinger picture wave function 

Consider the 3 + 1 dimensional quadratic theory, 

S = fdtd 3 x-[X 2 - (d t X) 2 -m 2 (t)X 2 }. (128) 

We work in a box of volume V = L 3 , with periodic boundary conditions. Define the modes, 
X n to satisfy the equation, 

X = £X n eW. (129) 

n 

The action, eq.(128), becomes, 

S = Jdt^[X n X. n + (^) 2 + m 2 (t)|X n | 2 ]. (130) 

Taking 

VVX n -► X n (131) 

gives the action for a single mode, 

S = Jdt[\X n \ 2 -u 2 (t)\X n \ 2 }, (132) 

with a time dependent frequency 

UJ \ t )=w 2 (t) = (^) 2 + m 2 {t). (133) 

Next, we calculate the wave function for the ground state, eq.(26), in the X, and X de- 
scriptions. 



20 This includes the transformation of the spin connection. 
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In the position space representation the ground state wave function is given by, 

rj}(x,i)=<x,t\0>, (134) 

where \x,t > is an eigenstate of the operator X(t), eq.(21). By definition, \x,t > satisfies the 
condition, 

X\x,t >— x\x,t > . (135) 

In this representation, P the canonically conjugate variable to X, is the operator, 

P = -id x . (136) 

Now from the definition of the ground state, eq.(26), and the expression for X, P, in terms 
of the creation and annihilation operators, eq.(21), eq.(22), it follows that the ground state 
satisfies the condition, 

P|0 >= (j)*X\0 > . (137) 

From the properties discussed above it then follows that the wave function tp(x,t), eq.(134), 
satisfies the equation, 

-id x ^{x,t) = {^Yxi){x,t). (138) 



This can be easily integrated to give, 

rl){x,t) = C{t)e^V"tK (139) 



Mr°t] 



The time dependent function C(t) is determined by requiring that Schrodinger's equation, 

~& x il>(x, t) + ^ 2 (t)^(x, t) = td^(x, t), (140) 

is met. It is straightforward to see that this gives, eq.(27). 
The function, /(£), is defined in eq.(16). At small t, 

H^i-coot) ~ iN u (-oj t) - c 1 (u t)- v , (141) 

where C\ is a constant. This leads to eq.(29). 

The probability density to find the system between x and x + dx is given by \ip(x, t)\ 2 . From 
eq.(27) this takes the form, 

\rl>(x,t)\ a = L±e- [ -W- ] . (142) 

We have used the fact that f(t) solves eq.(19), and has the asymptotic value, eq.(20). This 
means that the Wronskian, which is time independent, is given by, 

/(*)/•(*) " t(t)f(t) = 2iu . (143) 
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We see from eq.(142) that the probability density is a Gaussian with a width, J (Ax) 2 ~ |/|. 
This diverges as t — *■ 0, since |/| blows up, eq.(29). Thus the probability density gets more and 
more uniformly spread out as one approaches the singularity. 

The phase factor, eq.(30), arises from the limiting form of f(t) given in eq.(29). 

Next we turn to the wave function for the ground state in the X variable. The steps are 
analogous to those above. The X and X variables are related by, eq.(8). From the Lagrangian, 
eq.(14), it follows that the conjugate momentum, P, is given by, 

P = e~*X = e- $/2 (P + |x), (144) 

where we have used the relation, X = P. 
The relation, eq.(137), then leads to, 

j p|0>= e -*((4* + |)X|0>. (145) 

Let \x,t > be eigenstates of X, satisfying the condition 21 . 

X\x, t >= x\x, t > . (146) 

The wave function in the \x > representation, ip(x,t), then satisfies the condition, 

-i9^ = e-*[(^)* + |]^. (147) 

In addition the Schrodinger equation which takes the form, 

* ' P + e —oo 2 ^ = id t i>, (148) 



2 dx 2 2 

must be satisfied. This leads to the solution, eq.(35). The probability density \ip\ 2 , eq.(36), is 
then obtained as in the discussion above leading upto, eq.(142). To understand the behavior 
of the phase factor discussed in eq.(37), we note that /(£) is defined in eq.(16). At small t it 
then follows from the behavior of the Neumann function iV,,, that, 

f(t) ~ Cl (-u t) l *-»(l + c 2 {-uo tf). (149) 



21 We also require that the completeness relation, 

is satisfied. A similar relation, with x, replaced by x, is also satisfied by the states, \x,t >. This tells us that 
the states \x, t > are related to the states, \x, t > introduced above, by the relation, 

|x,t>= e -*/ 4 | a; >= e -*/ 2 |x,i>. 
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Thus, 

f S) 1 1 - 97/4-Tl 

2c 2 u;^. (150) 



./,„ , $ l,l-2i/ + ^ , n _ 2j 



7' 2 t v 2 
Here we have used, eq.(ll). This leads to eq.(37), after noting, eq.(17). 

The expectation value for P can be calculated from the wave function eq.(35). Alterna- 
tively, for this purpose, we can directly work in the Heisenberg picture. From eq.(144), and the 
expression for X, P, in term of the creation and annihilation operators, eq.(21), eq.(22), we get 
that in the vacuum state, 

<p> =W /+ i /! ' (151) 

Using, eq.(150) and related discussion above, this leads to eq.(39). 

A similar analysis can be carried out for a coherent state, defined by, 

a\s >= a\s > . (152) 

This leads to a wave function, 

iH.t) = j e [>-'4[({r t |i e ^y^/^. (153) 

The extra terms, compared to the ground state wave function, which are dependent on a, are 
both well defined in the limit t — > 0. Thus this wave function has the same type of singularity 
as the ground state wave function. 

C Subleading Contributions to Energy 

In this appendix we calculate the subleading contributions to the energy. These contributions 
would be the dominant ones if < V > vanishes, as discussed in section 3.4, and need to be 
calculated to understand when the energy remains finite, as t — > 0. 

A subleading contribution arises from the kinetic energy term. From the wave function, 
eq.(43), we find that for p > 1 this is given by, 

< KE >~ e —G 2 < {V') 2 >~ (-£) (2 ~ p) , (154) 

and diverges if p > 2. For p < 1, since G(t) is small near t — 0, 

e ' - ' < I 1 2 



<KE>~— J dx\ip' \ 2 ~ (-t) p . (155) 



This does not diverge as £ — *• 0. 
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Another subleading correction arises due to a correction in the absolute magnitude of i/j 
which in turn leads to a correction in < V >. We write, 

if>(x,t) = e- iG(<)y(£ Vo(5)[l + Si(x,t)]- (156) 

Since we are interested in the corrections to the absolute value of ip we take Si to be real. From 
the Schrodinger equation we get, 

y["^(|) + 2G(t)V'Re{f) + G(t)V") = ^. (157) 

For p > 1 the second and third terms within the square brackets on the left hand side dominate, 
leading to, 



5x = [/rfte $ G(t)][^i?e(7) + — ] (158) 

Wo * 



*2^ D „yo^ , ^ 



t 2 [\/'i?e(^) + --]- ( 159 ) 



2(1 -p) L V VV 2 

We see that this goes like t 2 , as t — ► and does not diverge. 

For p <\ the first term on the left hand side of eq.(157) dominates, giving, 



5x = ±J ro A)[/<ft e »] (160) 

z Wo J 

1 ./ m (^)(_ t )d+p). (i6i) 



2(l+p) >o 

This term goes like (— £)( 1+p ) and also does not diverge. Since S\ vanishes as t — > 0, the resulting 
correction to < V > and therefore to the energy also vanishes. 

Thus the conclusion is that that except for the case where p > 2, in which case the kinetic 
energy itself gives a divergent contribution, it is enough to have < V > as defined in eq.(52) to 
vanish,to ensure that the expectation value of the energy stays finite. 

D Particle Production 

In this appendix we detail the calculation of particle production at the quadratic level in the 
case where e $ does not become zero at any point, but can become small. For this purpose we 
choose a dilaton profile of the following form 

e Ht) = g s \t\ p \t\ > e 

e *(*) = g s \ e \ p \t\ < e. (162) 

We will first perform the analysis for each individual momentum mode, Xj,. 
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As explained above, one should work in the variables X since these are the variables which 
have a finite limit as t — ► 0. The equation of motion for X k is 



d ( . 
~dt \ 






x k = o. 



(163) 



Clearly, it is convenient to work with a time variable r defined by 

dt dr 
As is standard, we will solve (163) separately in the regions t < — e 
then match X k and d T X k across t = ±e. With the profile given in (162), a solution which is 
purely positive frequency at t — ► — oo is given by 



(164) 
e < t < e and t > e and 



X k (t) 

x k (t) 



uJ t p+1 



A exp 



J eP(p+l] 



B exp 



t< -e 

co t p+1 



eP(p+l) 



e< £ < e 



X k (t) = {u> Q t) v [C H?\u t) + D H<?\u> t) 



t> e, 



(165) 



where v has been defined in (17). The Bogoliubov coefficients C and D will be determined by 
the matching conditions. After a standard calculation we get the following expressions for C 
and D : 



C 



Z7T 
T 



(u; e){cos 



/ 2uj e \ r 

vP + iy 



H^urfH&iuoe) + H^MH^oe 



(i) 



\P+l) 



sin 



r(2) 



ffiiMwbcjHii^wbc) - tf^W^W) } 



D 



in 



4 (c oe ){2cos (^j i^W^iM 



. /2cu e\ 
— sin 



(4L\(^)) 2 - (^W)) 2 }■ 



In deriving these we have used the following property for any Bessel function Z v (x) 

^-\x v Z v {x)\=x v Z v _ x (x). 
dx 

A straightforward calculation verifies the unitarity relation 



\C\ 2 -\D\ 



-1. 



(166) 



(167) 



168) 



Let us first consider the limit uj e <C 1. Using the standard expansions for the Hankel 
functions, 



2^sin(vrz/)r(l + z/) 



4(i/ + 1) 

x\ 2u T{l + u) ' 
2 T(l-u) 



1 



x 



4(1 - v) 



+ 0(x 4 ))}, (169) 
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we find, 

in 



c 



4 sin 2 (7n/) 2 2 ^ [r(l - v)]< 



-2 2 ^ ( "° e)1 ' P - 2^u cosfvr,) F(1 = "} + ^V 



1-1/ ~ K ' r(i + i/) p+i 

(170) 

Thus there is a qualitatively different behavior of the Bogoliubov coefficient for p > 1 and for 
p < 1 as (c^oc) - ► 0. When p > 1 the coefficients C and Z) both diverge in this limit (of course 
maintaining the unitarity relation). When p < 1 they both tend to finite limits 

Lim W(je ^o C = — icot(7Ti/) Lim Woe ^ -D = — ie~ %1TU cosec(iiv) . (171) 

This difference between the cases p < 1 and p > 1 is the Heisenberg picture manifestation of 
the behavior of the Schrodinger picture wavefunctional. 

The analysis performed above was with an abrupt modification of the dilaton profile. How- 
ever we expect that the u^t <C 1 behavior would continue to be similar for a smooth modifica- 
tion. 

In the above analysis there is a finite amount of particle production for every momentum 
mode for p < 1, independent of the value of u>o, in the limit uoe — > 0. It is interesting to 
estimate the total amount of energy produced. However, for this estimate we need to perform 
the calculation for a dilaton profile which tends to a constant at early and late times, in keeping 
with our overall scenario. For a smooth dilaton profile, the ultraviolet behavior (uoe ^> 1) is 
then expected to be exponentially damped, \C\ 2 ~ e _tJo<E , so that the total energy produced is 
finite. 

It should be emphasized again that all the considerations of this appendix relates to the 
quadratic approximation. As we have seen this is not a good approximation in our problem. An 
estimate of the total amount of energy produced in the real problem has to take into account the 
effects of interactions which become stronger at later times. This requires a lot more detailed 
knowledge of strong coupling physics. 

E The Milne Background 

In this appendix we analyze the behavior when the boundary theory lives in Milne space with 
metric, eq.(92), with a dilaton, 

e ,J> = ^|tanh(t)| v/5 . (172) 

The metric eq.(92), upto the overall conformal factor, e~ 2t , is, 

dv 
ds 2 = -dt 2 + + r 2 (dd 2 + sin 2 ddtf). (173) 

This is a space of constant negative curvature. Since the gauge theory is conformally invariant 
it is equivalent to consider it in the background metric, eq.(173), and with dilaton, eq.(172). 
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Below we first analyze a scalar field, with Lagrangian 
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S = Jd 4 xV^e-*{±(dX) 2 + ^RX 2 }, (174) 

in the background with metric, eq.(173), and dilaton, eq.(172). Thereafter turn to the gauge 
field. The analysis in the scalar field theory, eq.(174), is not identical to the gauge theory, but 
quite analogous. 

We can mode decompose the scalar field into modes which are eigenfunctions of the 3 dim. 
spatial Laplacian. From [12], we see that, for the metric, eq.(173), the modes satisfy the 
equation, 

V 3 y k = -{k 2 + l)y k . (175) 

The functions, y^, are normalized to satisfy the condition, 

/ d 3 xVhy k (x)yt,(x') = 5 3 (k - k'). (176) 

Here hij is the spatial part of the metric, eq.(173). We can expand the field, X, in these modes, 

X = J d 3 kX k (t)y k . (177) 

This gives rise to decoupled oscillators for each mode, with the Lagrangian, 

S= fdtd 3 ke-' s '{\X k \ 2 -k 2 \X k \ 2 }, (178) 



where we have used the fact that the Ricci scalar, R = —6, for the metric. eq.(173). We 
see that for each mode the Lagrangian, eq.(178), is essentially the same as eq.(31), with k 2 
being identified with u 2 . Since the dilaton asymptotically goes to a constant here, eq.(172), 
the Lagrangian for each mode reduces to that of a standard harmonic oscillator in the far past 
or future. 

Our discussion in section 5 then leads to the conclusion that the wave function has a phase 
factor which is singular as t — ► 0, in this case as well. The phase factor is given by, eq.(61). 
The potential energy, V[X], for the Lagrangian eq.(174) is, 

V [X] = J d 3 xVh[h ij ^d i Xd j X + — RX 2 ]. (179) 

For the case of the gauge field, in the background metric, eq.(173), with the dilaton, eq.(172), 
an analysis similar to that carried out here and in the previous appendix can be done. Once 
again choosing Coulomb gauge is convenient. In this gauge one finds that the wave function 
has a phase factor which near the singularity takes the form, eq.(71), with, 

V[A i ( y x)] = jfd 3 xV^hF ij F ij . (180) 

This phase factor diverges, leading to a singular wave function. 



22 The conclusions would be essentially the same without the curvature coupling term, j^RX 2 . 
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F Universal behavior near singularities 

In this appendix we discuss some aspects of the universality of Kasner like behavior near space- 
like singularities in the class of models we consider. We will consider ten dimensional metrics 
of the form (105) so that the 3+1 dimensional Ricci tensor R^ and the dilaton $(x) satisfies 
the equations (106). It is therefore sufficient to discuss 3+1 dimensional dilaton cosmologies. 

Consider an AdS cosmology where the 4-metric is a Bianchi-IX spacetime. Using a BKL- 
type argument [13, 16, 15] we show the spacetime near the singularity has Kasner-like behavior. 
Furthermore, the dilaton drives the system towards an attractor region, where all the exponents 
are positive pi > 0, through a finite number of Kasner oscillations. This analysis can be directly 
extended to all homogeneous spaces with the results that we either have no oscillations at all or 
the number of oscillation is finite. This suggests that the symmetric Kasner singularity is generic 
and independent of the spatial 3-geometry, being either flat or any of Bianchi homogeneous 
spaces. 

It is worth mentioning that with no dilaton, symmetric Kasner solutions do not exist and 
the canonical BKL analysis gives an oscillatory approach to the singularity, with transitions 
between distinct asymmetric Kasner regimes. 

Let us take the 3+1 dimensional boundary metric g^ u (x) in (105) to be of Bianchi-IX type, 
which has the following form 

ds\ = -dt 2 + (a 2 {t)l a lp + b 2 (t)m a mp + c 2 (t)n a rip) dx a dx (3 (181) 

where l,m,n are the three frame vectors e^e^e 3 , (for explicit form of the metric see, for 
example, [13] page 390). The spatial symmetry algebra here is SU(2). a,b, and c are three 
independent scale factors 23 . 

If we assume that the dilaton is spatially homogeneous, then <9 a $ = e"9 Q $ = 0, with do& 
nonvanishing. Decomposing the Ricci tensor along the frame, we then have, 

*K + H=>> 2 - < 183 > 



with constant curvature, Rij — \^ij ,R = 



If wc take equal scale factors a — b — c, the spatial metric becomes da = {dx{ + dx\ + dx\ + cos x dx dx ) , 

3 

2' 
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and the dilaton field equation is given by 

In h r\ 

= 0. (184) 




The above system of ODE's is very difficult to solve analytically but if we ignore curvature 
terms (i.e., terms in (182) with no time derivatives) one has a = t Pl ,b = t P2 ,c = t ps , e $ = t a , 

I> = 1, E^ = l-y, (185) 

as an approximate solution near the singularity at t — 0. 

In the usual BKL analysis, a — 0, which forces one of the pi to be negative. The negative 
Pi means that time evolution towards the singularity necessarily makes one of the curvature 
terms ( treated as a perturbation to the time-derivative terms ) dominate the others at some 
point, for instance a 4 ~ £~' Pl ' (if p\ < 0). This forces the metric to evolve and transit from 
one Kasner regime to another according to the following law 

"' "17^' " 2 "77^' "" "TTi^' (186) 

However, with a nontrivial dilaton, we can have one of the following situations: All pi > 0, 
in which case no transitions take place since curvature terms (perturbations) die off as we 
approach the singularity e.g. 

[(&2 _ c 2)2 _ fl 4] _ _ fl 4 _ t 4/3 ^ q 

near the singularity t — > 0. Other Ricci components have similar behavior. This means that 
the symmetric Kasner case with all pi — | is stable against these perturbations as we approach 
this dilaton-driven symmetric Kasner singularity and there is no forced transition to a distinct 
Kasner regime. The other possibility is that one of the PiS is negative, in this case, we can have 
a finite number of oscillatory transitions between different Kasner regimes. This occurs since 
with every transition a increases and as it reaches a specific value (a = 1) all the pi become 
positive (see (185), and then no further transitions occur. 

If we consider any of the other types of Bianchi spaces, we should replace the curvature terms 
in (182) with that of this Bianchi space. But these terms have no time- derivatives so they will 
not change the leading behavior of the solution. Furthermore, these curvature terms either die 
off as we approach the singularity, in which case we have no oscillation at all, or one of them 
(pi < 0) gets larger as we approach the singularity. This leads to a finite number of oscillations 
as in Bianchi-IX. Again the dilaton will drive the system to an attractor region where the 
oscillation stops. 
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We will now see that the dilaton in fact drives the system towards an "attractor" region 
given by Pi,P2,P3 > 0, through a finite number of oscillations. Once the system reaches this 
region, there is no further oscillation. 

The equations 

5> = 1, E^ = l-y, (187) 

i i " 

for a Kasner-like cosmology with dilaton e* = t a can be described by the following parametriza- 
tion 



1— x \J\ — a 2 + 2x — 3x 2 1 — x y/l — a 2 + 2x — 3a; 2 . 
Pi = x, p 2 = — ^— + - , P3 = — - , (1* 

in terms of pi,a. For a solution to exist, the radical being positive forces 



< Pi < • (189) 

The above range for p\ can be divided into three regions. In these regions the values of the p^s 
are permuted among each other, as in figure (1). To avoid redundancy one should constrain p\ 
to one region which we choose to be 

1 - V4 - 3a 2 2 - V4 - 3a 2 , . 

—„ < Pi < — a • 190 

3 o 

The square root here being positive implies a 2 < |, i.e. \a\ < -m ~ 1.1547. The lower limit 
on px becomes positive if 4 — 3a 2 < 1, i.e. a 2 > 1. At this point, x — 0, a 2 = 1, all pi > 0. This 
shrinking of the allowed space of {p{\ is a key difference from the case a = without dilaton. 

Now let us say the system starts with say p[ = rr* •' < 0. Then there is a transition to a 
new Kasner regime with p\ and a^ given by 



(1) 


-ri 0) 


(D ri 0) +2rf 0) 


rf> = 


d 0) + 2ri 0) 


1 - 


"l + 2ri 0) ' 


l + 2p< > 




" 1 + 2p[ 0) ■ 





(191) 

Now, p\ < means that a*- 1 -* > a/ -*, z.e. a increases under the Kasner transition. More 
generally, for p_ < and p + > being either of the other two positive exponents, this can be 
expressed as the iterative map 

(„+!) -p. („ + i) P+ +2p_ a " H92) 

l + 2p l L> J lH-2pL ] l + 2p™ 

for the bounce from the (n)-th to the (n + l)-th Kasner regime with exponents Pi,Pj. The 
fixed point of this transformation is a = 0, and it is unstable for p_ < (an iterative map 
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Figure 1: Here P3 and P2 are plotted as functions of pi for a = case. Notice how the regions 
— 1/3 < pi < 0, < pi < 2/3, 2/3 < p\ < 1 have the same values for p^s 



x n +i = f(x n ) has an unstable fixed point x* 
increase of a is small for small a, since 



f(x*) if f'(x*) > 1). Furthermore, the rate of 



«n+l 



dr. 



Or 



-2p_ 
l + 2p_ 



(193) 



Thus a system with near constant dilaton (a ~ 0) takes a long iteration time to "flow" towards 
the pi > attractor region (although the flow is ensured due to the unstable fixed point). For 
instance, with p\ = xq = 0.3, «o = 0.001, the system flows (initially slowly) to pi > after 15 
oscillations, with 0:15 = 1.0896. 

This flow towards the Pi > attractor region in {pj}-space can be seen geometrically: the 
intersection of the sphere J2 Pi — 1 — \ with the plane J2 Pi — 1 is a circle on the plane. Under 
the bounce iterations, the sphere radius shrinks and so the circle radius also shrinks until the 
circle lies entirely in the Pi > quadrant. 

The finiteness of the number of oscillations means that the bulk cosmology flows towards 
the pi > attractor region, driven by the dilaton. 

Before we close this section we would like to comment on the nature of the flow towards the 



attractor region. Inverting (192), we get p_ 



Jl 



(™+i) 



l+2p 



(n+l) 



and so on 24 . Thus we can trace back 



from the symmetric Kasner, giving (upto five iterations) the flow, 



1 9 33, 
'5' 35' 35 ' 



5 7 19, 
21' 21' 21' 



3 5 9. 

TT'TT'TT 



1 3 3, 

5' 5' 5' 



'I I b 
3'3'3' 



(194) 



» 



since \p_ | < i is always true, therefore p. 



(n+l) 



l+2p ( L " ) 



> 0. This means that for each of the other two 



distinct p n > 0, we can potentially trace back to a distinct p_ < 0. This gives a tree with two flows starting 
with a given point {p\ n }. Similarly at every previous point, the tree forks into two. 
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An alternative distinct flow to the same symmetric Kasner endpoint begins at ("if>H'i§)> 
merging with the above flow at (— ]-[,yx,yx)- With each step backwards, a decreases. This 
shows that there are multiple trajectories that get attracted to any of the points in the {pi > 0} 
attractor region, perhaps as for any attractor-like behavior. 

Furthermore we suspect that the flow exhibit chaotic behavior, i.e. small changes in the 
initial conditions give rise to drastic changes in the final endpoints. For example, consider 
changing the starting point for the flow (194) above by a small perturbation (by ^ ~ 0.014, 
i.e. a 7% change to the smallest exponent, — |). This gives: 

13 9 65 2 13 39 3 2 23 13 17 

^~70'35'7(r "*■ ( ^~lT'44'44' ) ^ ^~28' 7' 2*T ^ ^ 11 ' 22' 22' ' ^ ^ 

the flow endpoint being distinct from the symmetric Kasner. 

We have used rational Kasner exponents above for simplicity in illustration: more gen- 
erally, one expects that there exist 'nearby' (not necessarily rational) Kasner exponents p.i 
with nonconstant dilaton (a ^ 0) in the neighborhood of exponents with constant dilaton 
(a = 0). In this case, a small change in the exponents in the set {pi,a ^ 0} would give 
exponents in the set {pi,a = 0}, which latter set belong to the canonical BKL analysis and 
oscillate forever, thus exhibiting no attractor behavior. More explicitly consider exponents 
{Pi>P2,P3,ct = 0} corresponding to a non-dilatonic asymmetric Kasner cosmology which os- 
cillates indefinitely, and perturb infinitesimally as {p[ = p\ + e,p' 2 = P2,p' 3 = P3 — e}. Now, 
a 2 = 2(1 — J2p[ ) ~ 4(p 3 — pi)e 7^ 0, if p\ 7^ p 3 . This latter set {p' { , a^0} thus flows to the 
attractor region, while the former does not. These examples and arguments suggest that small 
perturbations to initial conditions apparently give rise to large departures from the endpoints, 
in other words, chaotic behavior. 
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